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ABSTRACT 



Context. Disk galaxies are the most frequent objects among the galaxy population in the local universe. However, the formation of 
their disks and substructures - in particular their bars - is still a matter of debate. 
' r"| I Aims. We present a physical model of the formation of J033239. 72-275 154.7, a galaxy observed at z = 0.41 and characterized by 

O ^, a big young bar of size 6 kpc. The study of this system is particularly interesting for understanding the connection between mergers 

I . and bars as well as the properties and fate of this system as it relates to disk galaxy formation. 

P ■ Methods. We compare the morphological and kinematic properties of J033239.72-275 154.7, the latter obtained by the GIRAFFE 

spectrograph, to those derived from the merger of two spiral galaxies described by idealized N-body simulations including a star 
C/) , formation prescription. 

. Results. We found that the general morphological shape and most of the dynamical properties of the object can be well reproduced 

by a model in which the satellite is initially put in a retrograde orbit and the mass ratio of the system is 1:3. In such a scenario, a bar 
forms in the host galaxy after the first passage of the satellite where an important fraction of available gas is consumed in an induced 
burst. In its later evolution, however, we find that J033239. 72-275 154.7, whose major progenitor was an Sab galaxy, will probably 
become a SO galaxy. This is mainly due to the violent relaxation and the angular momentum loss experienced by the host galaxy 
during the merger process, which is caused by the adopted orbital parameters. This result suggests that the building of the Hubble 
sequence is significantly influenced by the last major collision. In the present case, the merger leads to a severe damage of the disk of 
the progenitor, leading to an evolution towards a more bulge dominated galaxy. 
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■ 1, Introduction of the fraction of actively star forming galaxies (including lu- 
^ ■ minous IR galaxies, LIRGs), of the fraction of galaxies with 
k> The formation of disk galaxies remains an outstanding puzzle peculiar galaxies (including those with compact morphologies) 
/S in contemporary astrophysics (see Mayer et al. 2008 for a re- etc... This has been followed by simulations of gas rich mergers 
b ' view). According to hierarchical models of structure formation, (Springel & Hernquist 2005; Robertson et al. 2006; Governato 
■ - - mergers and interaction of galaxies are an essential ingredient et al. 2007; Hopkins et al. 2008), evidencing that under certain 
of galaxy formation and evolution. Earlier works and numerical conditions, a disk may be re-built after a merger. Indeed such a 
simulations show that the remnants of mergings of purely stel- conclusion had been already reached by Barnes (2002), but these 
lar progenitors are more likely to be elliptical galaxies (Toomre simulations are less convincing since they do not include a pre- 
1977; Barnes 1988; Barnes & Hernquist 1992; Hernquist 1992; scription for star formation. Lotz et al. (2006) have also analysed 
Lima-Neto & Combes 1995; Balcells & Gonzalez 1998; Naab a large suite of simulated equal mass rich mergers and find that 
et al. 1999 etc.) and recent studies extended this result to gas- most merger remnants appear disc-like and dusty. Moreover, this 
rich progenitors (Springel & Hernquist 2005; Robertson et al. scenario is also consistent with results of other simulations and 
2006; Hopkins et al. 2008). Whether this is compatible, or in- semi-analytical models which claim that, without merger pro- 
compatible, with the fraction of disk galaxies present in the lo- eesses, most of galaxies and their host dark matter halos can- 
cal universe depends on the typical number of expected merg- not acquire the required angular momentum to form disks (see 
ers by galaxy (see for instance Kazantzidis et al. 2007). Such Peirani et al. 2004; Puech et al. 2007 and references therein), 
predictions seems to be incompatible with observations which The formation of bars is also a fundamental issue in the evolu- 
suggest that disk galaxies represent the majority (70%) of the tion of disk galaxies, particularly since it has now been shown 
galaxy population observed in the local universe (see Hammer that only about one third of them were in place at z = 0.8 (Sheth 
et al. 2005; Nakamura et al. 2004 and references therein). To et al. 2008). 

help resolve such issues, Hammer et al. (2005) suggested that Our aim here is to reproduce with numerical modelling the 

disks can be rebuilt during the encounters of gas rich spirals, general morphology (presence of a bar, of substructures, etc.), 

Indeed such proposition was guided by the remarkable coinci- the dynamical properties of the gas component (velocity fields), 

dence of the redshift increase, up to z=l, of the merger rate, the photometric properties of stars (colors, star formation rate, 
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etc..) of J033239.72-275 154.7, a galaxy located at z = 0.41 
for which we have high quality imaging and kinematics. This 
work will also provide useful input to disk formation models, 
since it gives information on potential progenitors of the present- 
day galaxy disks, as well as constraints on their formation (ini- 
tial orbital configuration, mass ratio of the system, etc..) and 
possible fate. This work is an off'spring of a VLT large pro- 
gram entitled IMAGES ("Intermediate-MAss Galaxy Evolution 
Sequence", Yang et al. 2008) which gathers high quality kine- 
matics for a representative sample of ~ 100 massive galaxies at 
z = 0.4 - 0.75 and with Mj{AB) < -20.3. Using the GIRAFFE 
spectrograph at the VLT, the kinematic properties of 65 of these 
galaxies, J033239.72-275 154.7 for instance, have been derived. 
This galaxy lies at z = 0.41, is classified as a merger from anal- 
ysis lead by Neichel et al. (2008) and presents a big, young bar. 
This bar has a size of 6 kpc, i.e. is quite big since 72% of barred 
galaxies from a sample of 2000 galaxies from the SDSS have a 
size lower than this (Barazza et al. 2008). This bar has also an 
extremely blue color, consistent with a starburst, i.e. with ages 
well below few 100s of Myr. This contrasts with many bars in 
the local universe which are known to include relatively old and 
red stars (Sheth et al. 2003). This paper is organized as follows: 
in Section 2 we present the general morphological and kinetic 
properties of J033239.72-275 154.7 and in Section 3 a short de- 
scription of the numerical modelling and of the results obtained. 
We conclude in Section 4. 

2. General properties of J033239.72-275154.7 

J033239.72-275 154.7 has been presented first by Yang et al. 
(2008). Table [T] summarizes its overall properties, including 
photometry, morphological parameters and kinematics measure- 
ments (provided by the IMAGES data basis and can be retrieved 
from Yang et al. 2008, Neichel et al. 2008 and Puech et al. 2008). 
This object has a stellar mass of 2.0 x lO^^M©, a K-band mag- 
nitude of Mk = -20.94 and shows a peculiar morphology and 
kinematics. Its center is dominated by an elongated structure, 
most likely a giant thin bar of semi-major axis 0.85 arcsec or 
about 6 kpc. Neichel et al. (2008) found that its color is roughly 
b-z=0.8, i.e. typical of starburst (see their figures 8 and 12). 
This bar is embedded into a diff'use region, which is probably 
a disk, as argued e.g. by features which look like tidal arms (one 
very blue at the upper right and another on the left). It is, nev- 
ertheless, very irregular, with several blue clumps. Roughly at 
its center, there is a strong light condensation which we will 
sometimes refer to as the core, and which is probably a small 
bulge Its color is roughly b-z=1.8, i.e. as found in galaxies of 
late type to Sbc (Neichel et al. 2008). The bar is asymmetric, 
extending spatially much more towards the upper-right side of 
the core. On the other side of the core, the bar is redder and 
a part of it disappears at UV wavelengths. The discrepancy in 
color between the two sides is roughly 0.3 magnitudes, as shown 
by the CDFS-GOODS observations (b band observations, thus 
UV rest- frame at z=0.4). On the bottom of the galaxy, there are 
two bright adjacent knots, which dominate the rest- frame UV 
light. These knots have the color of a pure starburst and could 
be hardly resolved. Due to the irregularity of the disk and to 
the existence of the bright tidal arms, it is difficult to measure 
accurately its position and inclination angle. From the shape of 
the outer isophotes we find that the photometric major axis is 
roughly in the upper-right lower-left direction and the inclination 
should lie between 28° and 35°. The star formation rate (SFR) 
is moderate (SFRir = 6.6, S FR2S00 = 4.2) and the dust does 
not aff'ect much the overall color. Note that SFRs have been esti- 
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Fig.l. B-V-z color map of J033239.72-275154.7 from 
HST/ACS (upper left panel) and distribution of the galaxy within 
the GIRAFFE grid (upper right panel). The lower panels show 
the velocity field map (lower left panel) and the associated cr- 
map (lower right panel). 



Table 1. Properties of J033239.72-275 154.7 



Multi wavelength photometry, stellar mass and SFR 

M28OO Wi Ml M.tellar SFR28OO SFR^ 

log(Mo) Mo/yr Mo/yr 

-19.56 -20.10 -20.62 -21.04 -20.94 10.31 4.2 6.6 

Morphology and Kinematics 

(from GIRAFFE measurements, Puech et al. 2008; Yang et al. 2008) 
V//flf Yhaif incldisk size of the bar 
km/s kpc deg. kpc 
30 3.5 35.4 ~6 



mated using the 2800A and 15 micron luminosities, respectively. 
Those have been converted to SFR using the method described 
in Kennicutt (1998). Luminosities have been estimated using in- 
terpolations between observed photometric points provided by 
HST/ACS, EIS and Spitzer/MIPS, using the method described in 
Hammer et al. (2001 & 2005). The ratio SFR of the SFR is con- 
sistent with Av=0.22 (standard extinction curve) which is very 
low. Finally, the electron density is modest, including at the very 
blue knots. 

The velocity field (VF) is obviously complex (see Fig. [T]). 
The kinematical major axis is almost parallel to the bar and it is 
off'set by more than one GIRAFFE pixel from the bulge, towards 
the prominent blue knots (in the bottom). The pixel scale is 0.52 
arcsec and the full width half maximum spectral resolution is 
23 km/s. Details on the GIRAFFE instruments can be found in 
Flores et al. (2006) and Yang et al. (2008). The small VF am- 
plitude is dominated by the two star-bursting clumps and by the 
blue arm/giant tidal tail. Even so, the velocity amplitude is rel- 
atively small, of the order of 30 km/s. The velocity dispersion 
map is nearly featureless at a value of 30-40 km/s over most of 
the galaxy. It has, nevertheless, a clear, localised sigma peak (50 
km/s), which is off'set by 2 GIRAFFE pixels from the bulge and 
could coincide (to within half a GIRAFFE pixel) with a small 
and relatively blue clump on the upper-right located roughly 
where the blue tidal arm joins the disk. The signal to noise ratio 
in this pixel is quite high (~ 50). 

All the evidence argues that strong interaction and/or merg- 
ing is at work in this system. Indeed, the dynamical axis is off'- 
centered and passes through the region where the two bright 
knots lie. Moreover, there are strong tidal arms, and the value 
of Vrot/cr is small (about 1-2 depending on the exact value of the 
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inclination). Indeed, the velocity dispersion is about 30-40 km/s 
over most of the galaxy. Given the VF values and the measured 
inclination, the rotational velocity is far below the value ex- 
pected from the Tully-Fisher (TF) relation, which is Vrot = 125 
km/s. This value has been taken from Hammer et al. (2007) who 
have carefully estimated the various samples used to derive this 
relationship. Given the spatial resolution, it is not easy to derive a 
rotation curve. Nevertheless, thorough modelling of the velocity 
fields have been done by Puech et al. (2008), who investigate the 
evolution of the Tully Fisher relation. Given its kinematical and 
morphological properties, it is thus unlikely that J033239.72- 
275154.7 can be a rotating disk hosting two giant HII regions 
because: (1) the two knots have a z band luminosity correspond- 
ing to one third of the total z luminosity (the relative photome- 
try has been done using apertures, scaled to the full width half 
maximum of the knot luminosity, using the IRAF/DAOPHOT/ 
package) ; (2) the dynamical axis is strongly off'-centered; (3) 
the value of cr is high and that of Vrot low. 

Could merging of the two bright knots with J033239.72- 
275154.7 be compatible with the formation of a giant bar with 
a relatively blue color? A/^-body simulations have shown that a 
small companion merging with a disk galaxy could, depend- 
ing on its density and its orbit, either destroy a pre-existing bar 
(Athanassoula 1999; Berentzen et al. 2003), or trigger its forma- 
tion (Walker et al. 1996; Berentzen et al. 2004). 

3. Simulations 

We briefly describe the numerical methodology used to model 
J033239.72-275 154.7. We use idealized N-body simulations of 
the merger of two spiral galaxies. One of the most difficult part 
of this work results from the huge number of free parameters in 
the initial conditions, for instance related to the mass ratio and 
the orbital configuration of the system or the star formation pre- 
scription. We have performed more than 100 simulations and we 
present in this section our fiducial modelling. First, as mentioned 
in section[2l the dynamical axis of the system is off'-centered and 
the velocity field strongly indicates that the satellite is falling 
down toward the main object. Then, by trying diff'erent orbital 
configurations, we found that retrograde orbits for the satellite 
are more plausible regarding the orientation of the tidal arms of 
the main galaxy (in the left and upper right parts) and the specific 
position of the accreted satellite. Second, diff'erent mass ratios of 
the system, diff'erent inclinations between the two orbit planes 
and diff'erent pericentric distances have been tested at the same 
time in order to match the amplitudes of the observed VF and cr 
map. For instance, since the VF amplitudes are small, the incli- 
nation between the two orbital planes should be small. It is also 
worth mentioning that we restricted our study to parabolic orbits 
(Khochfar & Burkert 2006) to reduce the number of free param- 
eters. Finally, diff'erent parameters related to the star formation 
recipes (see below) have also been tested in order to obtain re- 
sults consistent with the observations. 

3.1. Initial conditions and numerical method 

Our galaxies are composed of a spherical dark matter halo (with 
a Hernquist profile, Hernquist 1990), a disk, composed of stars 
and gas, and a bulge. In both objects, the disk and the bulge rep- 
resent respectively 15% and 5% of the total mass. The baryon 
fraction used here is slightly higher than the cosmic baryon 
fraction (~ 16.7%) derived by Komatsu et al. 2008 in order 
to compare to and match the results from Barnes (2002). It's 
worth mentioning that the gas fraction in the disk is 25% in the 



Table 2. Masses (M), gravitational softening lengths (e) and 
number of particles (N) of each component used in the simu- 
lations. 





DM 


gas 


star (disk) 


star (bulge) 


M,ost (lO^^Mo) 


6.60 


0.31 


0.93 


0.41 


M,at (lO^^Mo) 


2.20 


0.21 


0.21 


0.14 


6 (kpc) 


0.1 


0.2 


0.2 


0.1 


Nhost 


100000 


48000 


25714 


24000 


Nsat 


33333 


32000 


5714 


8000 



host galaxy, whereas it is 50% in the accreted satellite. These 
gas fractions are consistent with the estimation of Liang et al. 
(2006) Gavazzi et al. (2008) and Rodrigues et al. (2008). Note 
that the Gavazzi et al. (2008) measurements have been made 
in the local universe, while Liang et al. (2006) using the M- 
Z relationship, were able to estimate the evolution of the gas 
phase in z=0.6 galaxies, and found that, on average, galaxies 
at z=0.6 have two times more gas than at present days. The 
host galaxy has a stellar mass of 1.3 x lO^^M© and is more 
likely a Sab galaxy. Galaxies are created following Springel et 
al. (2005). Dark matter haloes have a concentration parameter of 
Chost = 14 and Csat = 15 for the host and the satellite, respec- 
tively, in good agreement with Dolag et al. (2004). Their spin 
parameter, defined by A = J \ E \^^^/GM^^^ where / is the an- 
gular momentum, E is the total energy of the halo and M is its 
mass, is Ahost = ^sat = 0-05 in order that the rotation curves of 
the galaxies are closed to the baryonic TF relation (McGaugh 
2005). Indeed, their maximum circular velocity is 140 km/s and 
98 km/s, respectively. In our fiducial model, the satellite is put 
on a retrograde parabolic orbit (Khochfar & Burkert 2006) with 
a pericentric distance Rp = 2.1 kpc and initial separation of 20 
kpc. The inclination between the two orbit planes is 15° and 
the galaxy spins are opposed. Note that the initial separation be- 
tween the two galaxies is quite small in order to reduce the high 
computational cost of all experiments. A higher distance would 
decrease the halos overlap at the beginning of the simulation, but 
this should not change the main results and conclusions of this 
paper. 

The simulation is performed using GADGET2 (Springel 
2005) with added prescriptions for cooling, star formation and 
feedback from Type la and II supernovae (SN). Approximately 
275,000 particles are used for the experiment and the masses 
(M), gravitational softening lengths (e) and number of particles 
(N) of each component involved are summarized in tabl^H 

The cooling and star formation (SF) recipes follow the pre- 
scriptions of Thomas & Couchman (1992) and of Katz et al. 
(1996), respectively. Gas particles with T > lO'^K cool at con- 
stant density for the duration of a timestep. Gas particles with 
T < 2 X 10"^ K, number density n > 0.1 cm"^, overdensity 
^Pgas > 100 and V.z; < form stars according to the stan- 
dard SFR prescription: dp^/dt = c^pgas/tdym where refers 
to the stellar density, tdyn is the dynamical timescale of the gas 
and is the SF efficiency. Assuming a constant dynamical time 
across the timestep, the fractional change in stellar density is 
^p*/p* = 1 - exp(-c*Ar/^dyn). For each gas particle, we draw 
a random number (r) between and 1 and convert it to a star if 
r < Apjp^. 

Instead of assuming 'instantaneous' energy injection, we in- 
clude the eff'ective lifetime of SN progenitors using the rate of 
energy injection Hsn- For this, we consider stellar lifetimes in 
the mass ranges 0.8 M© < m < 8.0 M© and 8.0 M© < m < 



4 



S. Peirani, F. Hammer, H. Flores, Y. Yang and E. Athanassoula: A giant bar induced by a merger event at z=0.4? 




Fig. 2. Time evolution of the projected star number density. The Hght blue arrows indicate the specific rotation of each galaxy, while 
the yellow dashed lines show the motion of the satellite. In panel 6, we superposed the GIRAFFE grid. Each frame is 40 kpc x 40 
kpc in size. 



80.0 Mq for Type la and Type II progenitors respectively. Using 
a Salpeter initial mass function for Type II SN gives: 



\0.24 



(1) 



where ^^a^ = 10^^ erg, is the mass of the stellar population 
and 3.53 < T < 29 Myr. For Type la SN, the heating is delayed, 
since they appear = 0.8 - 1 .0 Gyr after the onset of star forma- 
tion. Following de Freitas Pacheco (1998), the probability of one 
event in a timescale r after the onset of star formation is given 
by: 



erg.s 



(2) 



Eqns (1) and (2) are used to compute the energy released by 
SN derived from a star particle / (£"/). A fraction y of this energy 
is deposited in the j^^ neighbour gas particle by applying a ra- 
dial kick to its velocity with a magnitude Avy = ^(2wjyEi/mj), 
where wj is the weighting based on the smoothing kernel and mj 
is the mass of gas particle j. We note that all gas neighbours are 
located in a sphere of radius Rsn, centered on the SN progenitor, 
to avoid spurious injection of energy outside the SN's region of 
influence. In the following, we use the following standard val- 
ues: 7 = 0.1, 7?5A^ = 0.4 kpc and = 0.01. 

3.2. Results 

Figure (O shows the time evolution of the projected stellar num- 
ber density of the system. After 0.36 Gyr from the beginning of 
the simulation, the system has a stellar mass of ~ 1.92 x lO^^M© 
and presents, from a morphological point of view, a general 
shape and patterns similar to the observed ones. The disc of the 
more massive galaxy has not been destroyed yet, reflecting the 
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Fig. 3. The star formation rate obtained from our model (left 
panel) and projected distribution of newly formed stars at t = 
0.36 Gyr (right panel). The arrows in the left panel correspond 
to the SFR at the first encounter (1), at t = 0.36 Gyr (2) and at 
the final plunge (3). 



very early stage of the dynamical process, and is viewed in Figs[2] 
and[3]at an inclination of ~ 26°, as found from its outer isophote. 
The bar of the host galaxy has roughly the right orientation and 
length (~ 5 kpc) and the projected position of the satellite rem- 
nant is in rough agreement. However, many difficulties remain 
e.g. in reproducing the blue arm or tidal tail on the upper right 
and the fact that the companion should be split into two distinct 
bright knots. 

In the present scenario, the bar is forming after the first en- 
counter with the satellite, since at that time there is a very strong 
triggering. In good agreement with the observations, it is par- 
ticularly clear in the newly formed stars, since those, being on 
near-circular orbits, are more prone to the bar instability. Also in 
good agreement is the fact that the newly formed stars are also 
located in the satellite remnant (fig. [3]). Moreover, according to 
our model, J033239.72-275 154.7 is observed when the satellite 
is about to have the second encounter with the host. At this spe- 
cific time, the SFR is ~ 9MQ.yr~^ (fig-EJ) in good agreement with 
the observational estimation. 
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Fig. 4. Maps of the velocity field (on the left) and cr field (on the 
right) of the gas component obtained from our simulation. 
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Fig. 5. Mass profile decomposition along the major axis of the 
final product. The mass profile of the final product is show by the 
black line. The green dotted-dashed line, the blue dashed line 
and the red dotted line indicate the disk component, the bulge 
component and the sum of them, respectively. 

The VF and cr map of the gas component derived from our 
numerical model are represented in the GIRAFFE format in 
Fig HI The determination of the radial velocities has been real- 
ized at the orientation of Figs [2l and [3l The observed VF and 
its amplitude are faithfully reproduced and show that the initial 
inclination between the two orbital planes is consistent. The cr 
map is relatively flat around a mean value of 35 - 40 km/s, in 
good agreement with the observations. It presents two maxima 
of roughly the correct value, but located at positions other than 
those observed. However, these two maxima are located in re- 
gions which are not well resolved in gas particles and thus the 
resulting pixel values are not significant. 

After 4.26 Gyrs, corresponding to the present time z = 0, the 
final object is characterized by a bulge-to-disk ratio close to 1 .04, 
as derived from the mass profile decomposition along the major 
axis of the final product (see Fig. [5]). For this decomposition, we 
used standard models for the disk and bulge, namely an exponen- 
tial disk (Sersic profile with index n = 1) and a bulge following 
de Vaucouleurs law (Sersic profile with index n = 4), in which 
the free parameters are the bulge and disk radius and the bulge 
and disk flux. The rotational support parameter Vmaj/o'c = 1-62, 
where Vmaj is the maximum of the major axis rotation speed 
and (Tc the central velocity dispersion. This latter is estimated by 
considering all star particles in a sphere of radius 0.5 kpc. These 
previous values obtained highly suggest the formation of a SO 
galaxy. 

4. Discussion and Conclusions 

In this paper, we presented photometric and kinematic data of 
J033239. 72-275 154.7 and made a detailed comparison with an 



A/^-body simulation. A first interesting and important result is that 
among all the simulations that we have realized with specific or- 
bital parameters and progenitor mass ratios, the combination of 
parameters such as the star formation history, the VF and cr map 
put strong constraints on models and can remove the degeneracy 
of good candidates based on morphological criteria only. 

We found that our fiducial model can roughly reproduce the 
general morphological shape and the total stellar mass of the 
object as well as the observed SF rate, VF and intensity of the 
cr map. This specifically includes the reproduction of the giant 
young bar, its location and shape, the relative projected location 
of the companion, the overall morphological shape of the galaxy, 
its low rotation and its off'-center dynamical axis. The model 
predicts a slow-down of the velocity of the host galaxy due to 
exchanges of angular momentum with the companion (see Fig. 
[2]). However, using such idealized simulations, we were unable 
to reproduce many features, such as the blue arm or tidal tail 
on the upper right part, the morphological patterns of the two 
bright knots of the satellite remnant and the location of the cr 
peaks. Such discrepancies with observations can be easily un- 
derstood as due to the huge parameter space to be covered by 
the simulations. These include the properties of the progenitors 
(halo :disk:bulge: gas mass ratios and relative extents, shape of 
their density profiles and kinematics), the geometry of the en- 
counter and of the viewing, the gas properties, the SF, cooling 
and feedback modeling. On the other hand, the discrepancies 
with observations may also have astrophysical origin and/or be 
due to a more complex formation scenario. Several such possi- 
bilities - such as multiple encounters or progenitors with specific 
properties - come to mind, but would increase substantially the 
already too large free parameter space. 

In spite of these difficulties, it's encouraging to note that our 
simple numerical modelling is able to build a consistent picture 
of the formation of J033239. 72-275 154.7. This system is consis- 
tent with having been formed from a merger of two objects with 
a mass ratio 1:3. The simulation indicates that a bar is forming 
in the host galaxy after the first passage of the satellite where an 
important fraction of available gas is consumed in the induced 
burst giving a plausible explanation of the observed blue colors 
of the bar and of the satellite remnant. Moreover, both VF and 
cr map derived from the simulation match with the observational 
values and thus support this scenario. In its later evolution, we 
found that J033239. 72-275 154.7 may become a SO galaxy, as 
suggested by the results of our simulations. This is mainly ex- 
plained by the fact that the host galaxy experiences a violent re- 
laxation and looses some angular momentum during the merger 
process due to the retrograde orbit of the satellite. Moreover, by 
loosing some angular momentum, some of the gas can sink to- 
ward the center of the galaxy where it can be converted into new 
stars and then accelerate the growth of the bulge. 

To finish, it is interesting to point out that recent numeri- 
cal studies have demonstrated that gas-rich mergers can produce 
remnant disks provided strong feedback processes and also if 
both stellar and gas component do not experience a significant 
angular momentum loss (Springel & Hernquist 2005; Robertson 
et al. 2006; Governato et al. 2007, Hopkins et al. 2008). We pro- 
pose in this work an extreme example where i) the initial orbital 
configuration of the merger event do not permit to satisfy the 
latter criteria and ii) the progenitors of the system (in our nu- 
merical model) may not be sufficiently gas rich to reform a sig- 
nificant disk according those past studies. However, the present 
results are consistent with the previous ones in the emphasis on 
the fundamental role played by the last major event in building 
the Hubble sequence (see e.g. Hammer et al., 2005 and 2007). 
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It would thus be interesting to determine the frequency of each 
orbital configuration during a merger event and to compare with 
the galaxy population in the local universe to confirm this galaxy 
formation scenario. 
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